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ABSTRACT 

The physical structure of a nuclear flame is a basic ingredient of the theory of Type la super- 
novae (SNIa). Assuming an exponential density reduction with several characteristic times we have 
followed the evolution of a planar nuclear flame in an expanding background from an initial density 
6.6 x 10 7 g.cirT 3 down to 2 x 10 6 g.cirr 3 . The total amount of synthesized intermediate-mass elements 
(IME), from silicon to calcium, was monitored during the calculation. We have made use of the computed 
mass fractions, Xjme, of these elements to give an estimation of the total amount of IME synthesized 
during the deflagration of a massive white dwarf. Using X IME and adopting the usual hypothesis that tur- 
bulence decouples the effective burning velocity from the laminar flame speed, so that the relevant flame 
speed is actually the turbulent speed on the integral length-scale, we have built a simple geometrical ap- 
proach to model the region where IME are thought to be produced. It turns out that a healthy production of 
IME involves the combination of not too short expansion times, r c > 0.2 s, and high turbulent intensities. 
According to our results it could be difficult to produce much more than 0.2 M of intermediate-mass el- 
ements within the deflagrative paradigma. The calculations also suggest that the mass of IME scales with 
the mass of Fe-peak elements, making it difficult to conciliate energetic explosions with low ejected nickel 
masses, as in the well observed SN1991bg or in SN1998de. Thus a large production of Si-peak elements, 
especially in combination with a low or a moderate production of iron, could be better addressed by either 
the delayed detonation route in standard Chandrasekhar-mass models or, perhaps, by the off-center helium 
detonation in the sub Chandrasekhar-mass scenario. 

Subject headings: supernova: general, explosive nucleosynthesis 



1. Introduction 

A distinctive feature of SNIa spectra near maxi- 
mum light is the presence of intermediate-mass ele- 
ments from Si to Ca moving at velocities 10000-16000 
km.s -1 . In particular the presence or not of the Sill 
absorption feature at 6150A is basic to identify a 
Type I supernovae as SNIa. According to observa- 
tions about 0.2-0.4 M of IME have to be synthesized 
during the explosion. An amount of IME larger than 
O.2M is also needed to explain the observed abun- 
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dances of these elements in the solar neighborhood 
(Iwamoto et al. 1999). Therefore, any attempt to build 
a satisfactory theoretical model of these explosions 
has to give the right amount of intermediate-mass ele- 
ments. 

In this respect the underproduction of intermediate- 
mass elements found in current pure deflagration mod- 
els of Type la supernovae is a long standing problem of 
theoretical modelization of these explosions. On fun- 
damental basis we know that these elements have to 
be synthesized below densities p 5 x 10 7 g.cirT 3 in 
an expanding medium. Nevertheless the resolution 
achieved by the best multidimensional models carried 
out to date is far from adequate to solve unambigously 
this problem. In this work we have taken another route 
to address this issue. Instead of performing a multi- 
dimensional calculation incorporating the unresolved 
scales as a subgrid we have simulated the evolution 
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of an expanding nuclear flame in planar aproxima- 
tion with very good resolution. We have then used 
the results to put limits to the amount of IME syn- 
thesized in a deflagration. The usefulness of our ap- 
proximation requires that the synthesis of elements 
within the mass range 28 <A< 40 takes place, for 
the most part, during the fiamelet regime. In such a 
regime the flame, although highly wrinkled and dis- 
torted, still preserves its identity. Unfortunately the 
precise value of the minimum density at which the 
fiamelet hypothesis still holds is not known. There 
are indications that transition from the fiamelet to the 
distributive regime in the combustion diagram could 
begin around 2 x 10 7 g.cirT 3 (Niemeyer et al. 1999; 
Zingale et al. 2005), and we can take it as a reference. 

Keeping that constraint in mind we have proceeded 
in the following way. First, the detailed structure and 
evolution of a nuclear flame propagating through an 
expanding background has been solved for a reason- 
able set of characteristic expansion times and the nu- 
cleosynthetic yields calculated using a simplified a- 
network of 14 nuclei, § 2. With that information we 
provide in § 3 a rough calculation of the total amount 
of IME synthesized in a typical deflagration of a white 
dwarf as a function of the effective surface of heat ex- 
change between ashes and fuel at the densities of inter- 
est, the turbulence intensity and the characteristic ex- 
pansion time. According to our results a subsonic de- 
flagration is able to produce more than 0.2M Q of IME 
although only in a limited region of parameter space. 
This point is discussed in the chapter devoted to the 
conclusions, § 4, as well as possibles issues and limi- 
tations of our work. 

2. Production of Intermediate-mass elements in 
an expanding laminar flame 

From the nucleosynthetic point of view IME are 
the main product of incomplete Si burning at Tg < 
5 followed by a-rich freeze-out until the quench- 
ing of nuclear reactions. During the explosion of a 
massive white dwarf such combustion regime is at- 
tained when the density ahead the flame drops below 
- 5 x 10 7 g.cirT 3 in a medium which is rapidly ex- 
panding. As density decreases, the laminar velocity of 
the flame goes down, its width enlarges and the nu- 
cleosynthetic products are enriched first in IME and 
later in 24 Mg and 16 0. Finally the flame becomes a 
thermal wave and no more nuclear combustion takes 
place: the flame dies. Even though the thickness of the 



zone where IME are synthesized in the planar case is 
already really narrow, around half kilometer (see be- 
low), the surface of the flame in a real situation is so 
distorted by hydrodynamical instabilities that the ef- 
fective volume might be enough to account for the ob- 
served amount of IME. One important consequence is 
that, once developed, the hydrodynamical instabilities 
make the effective burning velocity independent of the 
laminar velocity of the flame (§ 3). 

We have devised a numerical experiment to inves- 
tigate whether or not the adequate amount of IME can 
be produced by the above mechanism. First of all 
we have made an estimation of the total amount of 
chemical species generated by a planar flame moving 
through a rapidly expanding media. Although in the 
end the effective turbulent burning velocity will be- 
come independent of the laminar value such calcula- 
tion is worthwhile because: 1) it provides the mass 
fraction of the isotopes of interest to be used in the 
next section when turbulence effects are taken into ac- 
count, 2) It shows the futility of trying the synthesis 
of IME in a real explosion when the flame only moves 
at its laminar speed, and 3) knowing the physics of an 
expanding planar nuclear flame is an interesting topic 
by itself as it represents a natural extension of the work 
by Timmes & Woosley (1992) (TW afterwards). The 
ouptut of the expanding planar flame calculation gives 
a fundamental quantity: Aitijme, the mass of IME syn- 
thesized in a column with 1 cm 2 section. For a given 
amount of burnt material Am/,„„ !( we compute the mass 
fraction of IME, Xjme synthesized below a fixed start- 
ing density po for different expansion rates. Later on (§ 
3) we make use of Xime to get an estimation of the to- 
tal amount of IME ejected during a typical Type la su- 
pernova event. The quantity Am IM E was calculated by 
following the propagation of the nuclear flame in pla- 
nar geometry as explained in Cabezon et al. (2004). 
Basically the integration scheme uses an a-network 
with 14 nuclei and implicit thermal coupling with the 
energy equation to follow the combustion under iso- 
baric conditions. The resulting method is robust and 
stable, allowing to properly handle the rear tail of the 
flame where most of the elements of interest were syn- 
thesized. In addition an exponential reduction in den- 
sity ahead the flame, with a characteristic time t c , is 
imposed. The resulting set of equations are: 

^ = x r « Y k Y ' - Yi n J YiY J + Z AmYm - AiYi (1) 
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where r t j stands for particle reaction rates, A t is for 
photodisintegrations and the rest of the symbols have 
their usual meaning. The thermal conductivity cr was 
taken from Appendix 1 in Khokhlov et al. (1997) and 
we have taken an uniform initial density of the fuel 
po = 6.6 x 10 7 g.cnT 3 . The precise value of po is irrel- 
evant for the results as long as it is larger than the den- 
sity at which IME begin to be produced. At each step 
we enforce the temperature of the unburnt material, lo- 
cated far enough from the flame front, to evolve adia- 
batically. The instantaneous value of T and p far ahead 
the flame sets the value of P exl in Equation (4). Across 
the flame the evolution is governed by heat difussion 
and binary 12 C + 12 C reaction, whereas in the burned 
zone far from the flame front the incinerated material 
relaxes to either Fe-peak elements or to IME if den- 
sity becomes low enough. The integration scheme is 
a standard lagrangian code with enough mass-shells 
to resolve the structure of the flame. At the begin- 
ning of the calculation a small group of shells were 
put into nuclear statistical equilibrium at constant pres- 
sure. Very soon a steady nuclear flame developed and 
propagated through the cold fuel (composed by equal 
parts of carbon and oxygen). As the flame progressed 
we kept a constant number of mass-shells ahead the 
flame (around 100) in order to save computer time. Fi- 
nally the calculation was stopped when the tempera- 
ture in the burned zone dropped below two billion de- 
grees and all nuclear reactions virtually quenched. At 
the end of the calculation there were typically between 
2 x 10 4 to 10 5 shells in combustion, depending on the 
characteristic expansion time t c . It was necessary to 
calculate several millions of time steps in order to fol- 
low the whole process. 

In Figure 1 it is shown a snapshot of the mass frac- 
tion profile at t=0.1523 s, which is the time at which 
total IME abundance freezes in our reference model 
with t c =0.1 s. As we can see 28 Si as well as other 
IME with higher atomic mass were synthesized in a 
tiny layer, lesser than 0.3 km, which is more than 
one order of magnitude below the resolution of the 
best three-dimensional simulations carried out so far. 



Moreover because during the explosion the flame is 
shattered by turbulence that length represents an upper 
limit of the size of pockets of ashes made of elements 
belonging to the silicon family. Elements lighter than 
silicon such as magnesium and neon were produced at 
lower density in a very narrow region behind the flame 
edge. In Figure 2 we depict the evolution of the total 
synthesized mass of elements between 28 Si to 40 Ca for 
t c = 0.05, 0.1, 0.2 and 0.3 s as a function of the fuel 
density. From these figures it is clear that for plausible 
expansion times the IME are synthesized in the den- 
sity range 2 - 5 x 10 7 g.crrT 3 , as it was recognized 
many years ago. In particular, we see that 95% of 
IME were already synthesized at densities higher than 
2 x 10 7 g.cnT 3 . In fact, more than a half of their final 
mass was produced when the density of the fuel was 
still above 4 x 10 7 g.cnT 3 . This is an interesting fea- 
ture because it suggests that important nuclei for SNIa 
spectroscopy, such as silicon and calcium, are prob- 
ably synthesized during the so-called flamelet regime 
of turbulent combustion, much more easier to handle 
than the distributive regime. Below 2.5 x 10 7 g.cnT 3 we 
found an increase of direct carbon and oxygen com- 
bustion products, neon and magnesium in our network, 
in a region of size ^ 0.05 km. However in a real phys- 
ical situation the final abundance of elements below 
silicon might be altered if combustion enters the dis- 
tributive regime. The velocities of the flame, fitted us- 
ing the total amount of burnt fuel (see Eq.7 below), are 
in reasonable agreement with those obtained by TW 
at different densities. We obtain slightly lower values 
because the limited nuclear network used to carry out 
the calculations. For example, we obtained v=0.14, 
0.23, and 0.33 km.s -1 at p-i-2, 3 and 4 respectively, 
to be compared to those calculated by TW, V7-iy=0.19, 
0.32, 0.46 km.s -1 at the same densities. The instanta- 
neous width of the flame was also consistent with TW 
calculation. The detailed final yields of the 14 nuclei 
for several characteristic expansion times are shown 
in Table 1. As we can see, there is a monotonic in- 
crease of the final mass of the most abundant species 
with the expansion time t c . For t c = 0.05 s we got 
2.95 x 10 10 g.cnT 2 of silicon whereas for t c = 0.3 s 
the synthesized mass was 1.16 X 10 11 g.cnT 2 . The 
same trend is true for the size of the region where 
these elements are significatively produced: ^0.15 km 
(t c = 0.05 s) to =s 0.75 km (t c = 0.3 s). Even in the 
last case the width of the IME layer is much smaller 
than the size of the burned region, typically 2000 km. 
Thus in the lagrangian one-dimensional picture it can 
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be regarded as a very thin layer on top of a volume 
made basically of Fe-peak elements. 

The mass of fuel consumed during the propagation 
of an expanding planar flame starting at an initial den- 
sity po can be found analytically because the density 
dependence of the laminar velocity is known (TW): 



v = v — 



(5) 



From our results vo - 0.6 km.s -1 at p = 6.6 X 
10 7 g.crrT 3 , and n - 1.2. Using Equation 3 the mass 
of burned fuel per unit surface in an elapsed time Af is: 



Table 2 that can be also used to extrapolate AmiME be- 
yond 0.3 s, until t c =! 0.5 s if necessary, is: 



11 - -0 

Am /M£ = 6.5 x 10 t c 4 g.cm 



(8) 



As mentioned before, the mass of IME given by 
Equation 8 is independent of the particular value cho- 
sen for po in Equation 3 provided it is high enough 
to resolve the region where intermediate-mass ele- 
ments are synthesized. However their mass fraction, 
X/me, within the burnt material is a function of the 
starting density po in the exponential expansion. For 
po = 6.6 x 10 7 g.cm -3 it can be fitted using: 



kniburnt = PO VO 
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= 0.35 r r 



(9) 
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Taking Af — > oo in Equation 6 gives a good approxi- 
mation to the the total burnt mass per unit surface: 



Am burnt = p v 



(1+n) 



(7) 



Hence the amount of consumed fuel below the start- 
ing density po in the expanding background is simply 
proportional to the product of the laminar velocity and 
the characteristic expansion time. The laminar veloc- 
ity of the flame is rather well known within, at most, 
a factor two uncertainity. The value of the character- 
istic expansion time t c is probably dependent of the 
particular explosion dynamics. A rough lower limit 
for r e is the characteristic hydrodynamical time scale 
Th = 446/ y/p = 0.08 s at p = 3 x 10 7 g.cnr 3 which 
matches well also the behaviour of a particular de- 
flagrative one-dimensional model (Bravo et al. 1996) 
as shown in Figure 3. A reasonable range for t c is 
0.05s < t c < 0.3s once the flame has reached den- 
sities below 6 x 10 7 g.cm -3 , although higher values, 
nearing 0.5 seconds, can not be totally discarded. In 
Table 2 there is shown the total amount of burned mass 
as well as the masses of IME and nickel synthesized 
below po as a function of r c . As we can see the burnt 
mass scales linearly with r c , as suggested by Equation 
7. However, the amount of synthesized IME is not pro- 
portional to the characteristic expansion time because 
a longer value of t c leads to the formation of more iron 
as the nuclear system has more room to relax. A sim- 
ple fitting formulae to the values of Am /M£ given in 



Another remarkable feature shown in Table 2 is that 
the production of 56 Ni is no longer negligible, espe- 
cially for t c > 0.2 s. At the longest expansion time 
considered here, r c = 0.3 s, the amount of 56 Ni is prac- 
tically equal to that of IME. 

3. Consequences for Type la Supernova models 

A quantitative estimation of the total amount of 
IME ejected in a SNIa explosion may be done pro- 
vided the effective flame surface at densities p 
2-5 x 10 7 g.cirT 3 is somehow known. However, this is 
a difficult task that will have to be addressed in the fu- 
ture through high resolution multidimensional simula- 
tions. For the time being we give only a qualitative dis- 
cussion about this point taking advantage of the self- 
similar behaviour of the main physical agents respon- 
sible for the corrugation of the flame. There has been 
shown (Khokhlov 1995; Niemeyer & Woosley 1997), 
that a large increase in flame surface can only be at- 
tributed to two physical mechanisms: the Rayleigh- 
Taylor (RT) instability and turbulence. Fortunately, 
it turns out that both are self-similar phenomena that 
can be approximately described using three param- 
eters: the maximum and minimum scale-lengths, 
A max , A m i„ and the average fractal dimension D' of the 
corrugated flame surface. In this approximation the 
incinerated mass is given by: 



M bu mt = Inr 2 



D'-2 



Jo ' 



PV fl ame dt (10) 



where r is the average radius of the flame and 
v fiame stems for its laminar velocity. Even though 
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in the supernova the nuclear flame becomes highly 
distorted and shattered by turbulence and hydrody- 
namical instabilities the amout of fuel consumed per 
unit area given by the integral part of Equation 10 
should not be very different from Ami, urnt given by 
Equation 7 corresponding to the planar case ana- 
lyzed in the previous section. However many exper- 
iments in combustion chambers as well as numerical 
simulations dealing with turbulent combustion sug- 
gest that A m i„ is no longer an independent quantity 
(Peters & Franke 1989). Its local value self-adapts so 
that the effective burning rate is set by the integral 
length A max and the characteristic velocity, V, at the 
integral length-scale. Such interesting behaviour relies 
in the existence of an universal scaling law for velocity 
spectrum of the type: 

(£=)-(—)* <»> 

\ A m in j \ V fl ame J 

which for D'=7/3 becomes the popular Kolmogorov 
scaling-law for turbulence. Adopting the scaling rela- 
tionship given by Equation 11, both Mb urn t as well as 
the mass of the k-isotope, M* become independent of 
v flame- The ejected mass of the k-isotope can therefore 
be calculated using the value of computed in the 
previous section: 

M k =4ttF 2 V |JT 'pC)*) X k (12) 

where Atk is the elapsed time until the abundance of the 
respective nuclei freezes. For example, as stated in § 2 
the nuclear combustion is unable to synthesize nuclei 
beyond magnesium once the density has gone down 
=s 1.5 x 10 7 g.cirT 3 . Within our model of exponential 
density reduction this occurs once At IM E = 1-49 r c . 

Admittedly, the above discussion is a little artifi- 
cial, but necessary to develop the formalism, because 
finally neither Vfi ame nor D' appears in Equation 12. 
According to that equation the progression of the flame 
can be envisaged as a spherical flame front of radius 
f moving not at the laminar speed but at velocity V. As 
V is expected to be many orders of magnitude larger 
than v fiame the synthesized mass of nickel as well as 
the released nuclear energy are sufficient to account for 
the supernova light curve and energetics. Nevertheless 
a variation of the geometrical approach given by Equa- 
tion 12 could be more fruitful to address the particu- 
lar problem of the synthesis of IME (Woosley 2001). 
In this regard it might be more physical to interpret 



r as the average radius of the ashes, made basically of 
nuclear-statistical equilibrium elements (NSEE). The 
intermediate-mass elements are then produced in a 
tiny shell, almost a surface, on top of the huge vol- 
ume defined by the underlying ashes, over a length 
given by the product V At me- Unlike the synthesis of 
NSEE, which as bounded to a volume are reasonably 
given by Equation 12, surface effects are important to 
properly model the synthesis of IME. Assuming that 
the effective surface of NSEE behaves as a fractal in 
which those scale-lengths below the current thickness 
of the IME layer itself do not appreciably contribute to 
M/me, we can write: 




which either for D=2 or for (V AIime) — f reduces to 
Equation 12. Here, the fractal dimension D is intended 
to describe the geometrical properties of the NSEE ex- 
ternal surface in between the integral length-scale and 
the radius of the front, and has not to be the same as the 
dimension D' characterizing the front at lower length- 
scales. Finally, to compute Mime from Equation 13 it 
is more convenient to use the mass of the ashes made 
of NSEE, Mmsee, instead their average radius f. Then, 
solving for the integral, taking X IM e from Equation (9) 
and using At m E = 1 -49 t c the resulting expression is: 

M IME (M e ) = 0.253 | M "s*E(M ) j D/3 2A3 d y3-n t ii/4-d 

(14) 

where pi is the average density of the ashes in 
10 7 g-crrT 3 and Vs the average velocity of the effective 
flame velocity in 10 8 cm-s~'. 

The main results of a parametric exploration of 
Equation 14 are depicted in Figure 4. In this figure 
there is shown the synthesized mass of elements be- 
tween silicon to calcium as a function of the values 
adopted by the set {Vg, D, t c ). In computing Figure 
4 we have taken a fixed mass of ashes Mnsee(M q ) = 
0.8 and p-j — 7, a density which is slightly higher than 
the starting density during the exponential expansion. 
Such combination of M NS E e and p 7 leads to an average 
radius of the ashes f = 1.8 x 10 8 cm, in fair agreement 
to detailed hydrodynamical calculations. For instance, 
model W7 (Nomoto et al. 1984) gives approximately 
r 8 =s 2 at p-i = 4 similar as the one adopted in this 
work. Going to multidimensional models Ropke et al. 
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(2006) give an average radius of the burnt zone also 
close to ?s = 2 in their model c3_3d. 

First of all note that the case D=2 in Figure 4 al- 
most always leads to an underproduction of IME. Thus 
it becomes clear why those one-dimensional models 
which assumed a spherical burning front moving at the 
largest turbulent speed, V, did not succeed in synthe- 
sizing large amounts of IME (Niemeyer & Woosley 1997). 
Also, there is not a sufficient abundance of IME nei- 
ther for V 8 < 0.5 nor for t c < 0.05 s. The mass of 
these elements smoothly rises as D, V or t c increases, 
as expected. According to Equation 14 the most in- 
fluencial parameter affecting the nucleosynthesis of 
silicon family is the fractal dimension of the under- 
lying iron core, followed by the propagation velocity 
and the characteristic expansion time respectively. 

There are, however, some limitations which could 
prevent a large production of IME during the explo- 
sion. First, the fractal dimension defining the ashes's 
surface is thought to be closer to 2.3 than to 2.5 be- 
cause it is highly probable for turbulence to be of Kol- 
mogorov type. Actually, one can expect that the flame 
behaves as a multifractal, so the effective fractal di- 
mension would be somehow between 2.33 (character- 
istic of lower scale-lengths) and D=2.5 (characteristic 
of large scale-lenghts). On the other hand the local 
sound speed puts an upper limit to the largest turbu- 
lent velocity V. At the densities of interest the sound 
speed is 4 - 4.5 x 10 8 cm.s~' thus Vs -2 is almost a 
half of the sound speed. In that case the induction of a 
detonation because of the cumulative effect of pressure 
waves can not be ruled out (Plewa et al. 2005). Maybe 
more interesting are the constraints related to the dom- 
inant combustion regime: flamelet or distributed. Tak- 
ing A max = r = 1.8 X 10 8 cm, A m i„ ^0.1 cm (ap- 
proximately the flame thickness at these densities), 
v fiame - 5 x 10 4 cm.s -1 and D'=7/3 in Equation 1 1 the 
corresponding turbulent velocity is 6 x 10 7 cm.s~' . 
With that choice of parameters turbulent combustion 
will stay in the flamelet regime as long as the veloc- 
ity of the largest eddy is lower than 6 x 10 7 cm.s~'. 
Were this the case an inspection of Figure 4 reveals 
that, even in the most favourable case with t c = 0.3 s, 
the production of IME hardly arrives to O.2M . Con- 
sidering stronger turbulent intensities would move the 
more plausible region for producing IME just to the 
poorly known distributive regime. 

The influence of the characteristic expansion time 
t c is not as strong as that of Vs owing to the inverse 
dependence of X IME with t c (Eq. 9). Moreover, one- 



dimensional simulations indicate that the characteristic 
expansion time at the time combustion arrives to lay- 
ers with p =s 5 x 10 7 g.cirT 3 are probably not much 
larger than 0.2 s (see for example Fig. 3). Even though 
r c should be relatively model independent, owing the 
strong sensitivity of the degenerate matter against per- 
turbations, its precise value is set by the central den- 
sity at the moment of the explosion and by the nuclear 
combustion dynamics itself. 

Equation 14 suggests that Mime and M^see are not 



independent quantities, Mime k M 



d/3 

NSEE 



On the other 

hand some mass of Fe-peak elements are also pro- 
duced simultaneously to silicon during the expansion. 
Thus, the final quantity of Fe-peak, M Fe is slightly 
higher than M NSEE . Assuming that after freezing is 
X Fe =! 1 - Xjme in the IME layer it is easy to show that 
the total mass of the iron group is: 



M Fe (M e ) = M NSE e(M q )+M ime (M q ) 



' c 

0.35 



0.05 < T r < 0.3 



(15) 

A plot of Mime versus M Fe for several characteris- 
tic expansion times at constant D=7/3 and Vg=1.5 is 
given in Figure 5. As we can see the relationship be- 
tween Mime and M Fe is close to linear with increasing 
slope as t c rises. As a general rule we expect that de- 
flagrations with large luminosities (i.e a large mass of 
56 Ni) will also show a large production of IME. Nev- 
ertheless some attention has to be paid to the influence 
of the characteristic expansion time. For example, ac- 
cording to Figure 5, 0.2M o of IME can be synthesized 
in combination to 0.7M o of Fe-peak elements at t c = 
0.3 s or to 0.95M© at t c = 0. 1 s. Hence explosions with 
different maximum luminosities and global energetics 
could lead to similar amounts of IME. Moreover, it is 
reasonable to think that r c and M NSE e are inversely 
correlated: a large amount of M NS ee means a big deal 
of released nuclear energy thus leading to shorter char- 
acteristic expansion times. Such correlation suggest 
that some self -regulating mechanism could be at work 
reducing the dispersion in the production of IME. 

Finally, note that our assumption that ashes are 
enterily composed of NSEE has not to be necessar- 
ily true. Even though there is a current debate con- 
cerning this point, recent multidimensional simula- 
tions indicate that some unburnt carbon and oxygen 
could remain mixed with the incinerated material, 
(Garcfa-Senz & Bravo 2005; Ropke et al. 2006). In 
that case the correlation between the mass of IME and 
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Fe-peak elements will be different to that shown in 
Figure 5. 

4. Conclusions 

The amount of intermediate-mass elements ejected 
in a SNIa explosion is an important diagnosing tool to 
discriminate between a variety of theoretical models 
attempting to explain the nature of such explosions. 
One of the most popular models is the deflagration 
model, rooted in the subsonic propagation of a nuclear 
flame through the white dwarf interior. However, with 
very few exceptions, neither the one-dimensional sim- 
ulations nor modern three-dimensional calculations of 
this model have been able to give the right quantity of 
IME. Hydrodynamical calculations always stubbornly 
give a too low production. This fact is currently at- 
tributed to the low resolution achieved in multidimen- 
sional calculations. 

In the first part of this paper we have analyzed the 
synthesis of elements from helium to zinc during the 
propagation of a nuclear flame in an expanding back- 
ground. The simulation was carried out in the planar 
aproximation, thus the resolution was enough to re- 
solve the detailed structure of the flame. In this sense 
our work represents a continuation of that by Timmes 
& Woosley (1992). As a result we were able to get an 
idea of the size of the region where the IME were for 
the most part synthesized: between 0.15 km to 0.75 km 
depending on the characteristic expansion time. This 
width is hardly reachable by today multidimensional 
hydrocodes. It was also possible to give the amount 
of elements from silicon to calcium cooked behind 
the flame (Table 1). This information can be used by 
SNIa modelers to improve the nucleosynthesis output 
of their hydrodynamical simulations. An interesting 
result is that IME are massively produced in the range 
p = 2 x 10 7 - 5 x 10 7 g.cirT 3 . In fact, according to our 
calculations 80% of the IME were already produced 
before the density declined below 3 x 10 7 g.cnT 3 (see 
Fig. 2). This could have important consequences be- 
cause the usefulness of our results strongly relies in the 
basic hypothesis that the flame does not enter the dis- 
tributive regime above 2-3xl0 7 g.cirT 3 . Finally useful 
fitting formulae giving the mass of IME synthesized 
during the expansion as well as to their mass frac- 
tion, X/M£, as a function of the characteristic expan- 
sion time t c were also given (Eqs. 8,9). 

In § 3 we did use the information given in Table 
1 to estimate the quantity of IME ejected during the 



deflagration of a white dwarf taking advantage of the 
self-similar properties of the thermonuclear combus- 
tion in SNIa. We have taken a geometrical approach to 
address this problem in which the narrow region where 
the Si-peak nuclei are typically produced was treated 
as a tiny shell wrapping around a deformed surface on 
top of a volume made of Fe-peak elements. That ap- 
proach holds as long as the ashes form a well defined 
and not too dispersed pack limited by an irregular sur- 
face. Within our simple model the total amount of el- 
ements from silicon to calcium can be obtained from 
an analytical expression (Eq. 14) which uses four pa- 
rameters: the ratio between the mass and average den- 
sity of the NSE-ashes, the fractal dimension of their 
surface, the velocity of the turbulent combustion at the 
largest scales, and the characteristic expansion time. In 
consequence we have tried to delimitate the regions in 
parameter space able to give the amount of IME com- 
patible with observations (Fig. 4). 

As a general conclusion we state the difficulty 
of pure deflagration models to produce more than 
0.2M o of intermediate-mass elements. In this sense 
our phenomenological approach may help to under- 
stand why the underproduction of Si-peak elements 
has traditionally been a long standing problem of all 
sort of hydrodynamical simulations concerning Type 
la supernovae assuming the deflagrative paradigma. 
Still there remains a significant area in space param- 
eter showing a rich production of IME. In particu- 
lar the region defined by the unequalities t c > 0.2 s, 
D > 2.2 and V$, > 1 seems to be well suited for that 
purpose. The above lower limit in the fractal dimen- 
sion of the ahes's surface does not pose a problem 
as experiments dealing with turbulent combustion in- 
dicates a favoured value of D close to 2.3. As an 
aside, ID simulations that always assume D = 2.0 
on the resolved length-scales lead unavoidably to an 
underproduction of IME. Perhaps more demanding is 
the limitation imposed by the characteristic expansion 
time. Even though t c should not be too different from 
the hydrodynamical time scale, ^ 0.1 s at the densi- 
ties of interest, its particular value is model dependent. 
Hence expansion times with t c ^ 0.2 s are probably 
not unrealistic. Turbulent intensities characterized by 
Vs > 1 are also reachable during the explosion be- 
cause the rising velocity of the larger floating blobs 
of incinerated material is of the order, or even greater, 
than that velocity. Nevertheless, the turbulent cascade 
generated by these fast moving blobs could fragment 
the flame in very small pieces, shorter than the flame 
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thickness, invalidating our hypothesis that combus- 
tion always takes place in the flamelet regime. Thus, 
although our main goal was to highlight the favoured 
regions in parameter space to synthesize a good deal of 
IME, the presented results also underline the potential 
importance of the distributed regime for deflagration 
models, if they were to represent the main channel to 
normal Type la supernovae. 

Even if a large mass of IME were attained in de- 
flagration models of SNIa, our study still suggests 
another potential problem: the production of IME 
seems to increase monothonically with that of Fe- 
peak elements. Hence, it would be seemingly difficult 
for deflagrations to explain peculiar low-luminosity 
but energetic events like SN1991bg or SN1998de 
(Stritzinger et al. 2006). 
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AYA2002-04094 and also supported by DURSI of the 
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Fig. 1 . — Mass fraction profile of the most abundant 
species synthesized during the decay of a planar nu- 
clear flame due to expansion. A characteristic expan- 
sion time t c = 0.1s was assumed in this case. The 
elapsed time was t=0.1523 s, which corresponds to a 
fuel density of 1.45 x 10 7 g.cirT 3 . Below that density 
the total amount of IME practically freezes 
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Fig. 2. — Total mass per unit area of IME (from silicon 
to calcium) synthesized during the isobaric propaga- 
tion of the flame as a function of the instantaneous den- 
sity value in the cold fuel for four different characteris- 
tic expansion times. At p = 2 x 10 7 g.crrT 3 practically 
all IME have already been produced. 
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Fig. 3. — Evolution of density, radius and lagrangian- 
mass coordinate (in 10 7 , 10 8 and 10 33 c.g.s units re- 
spectively) of the one-dimensional model used to com- 
pute some of the inputs of Figure 4. The figure de- 
picts the value of these variables just at the moment 
the temperature of the corresponding mass-shell be- 
comes higher than one billion degrees. Thus, crosses 
represent the density of different mass-shells as they 
are reached by the flame. The solid line is an exponen- 
tial fit to the numerical data using a characteristic time 
t c = 0.095 s. We have checked that the maximum tem- 
perature achieved by these shells is in the right range 
to give intermediate-mass elements. 
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Fig. 4. — Total amount (in M ) of elements from sil- 
icon to calcium synthesized during the expansion of a 
massive white dwarf due to the propagation of a turbu- 
lent nuclear flame. The region where these nuclei are 
produced is approximated by a fractal of dimension D 
and average radius f = 1.8 x 10 8 cm. The quantity 
of IME is then given as a function of D, the effective 
propagation velocity of the flame V e /f, and the charac- 
teristic expansion time r e . 
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Fig. 5. — Mass of IME as a function of the ejected 
mass of Fe-peak elements for several expansion times. 
A constant fractal dimension of the iron core surface, 
D=7/3, was assumed as well as a constant turbulence 
intensity, V= 1.5 X 10 8 cm.s -1 . 
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Table 1: Final masses (g.cirT 2 ) of a-elements produced 
in a laminar expanding nuclear flame for 
several characteristic expansion times 



Nuclei 


t c = 0.05 s 


T e = 0.1 s 


t c = 0.2 s 


t c = 0.3 s 






-mass (g.cirT 2 )- 






4 He 


4.1 10~ 2 


1.62 10" 1 


4.25 10 5 


2.7 10 7 


12 C 


3.3 10 8 


3.48 10 8 


3.72 10 8 


3.9 10 8 


16 Q 


4.9 10 9 


7.7 10 9 


1.23 10 10 


1.6 10 10 


20 Ne 


3.5 10 7 


4.4 10 7 


5.60 10 7 


6.9 10 7 


24 Mg 


1.9 10 9 


3.2 10 9 


5.12 10 9 


6.9 10 9 


28 Si 


2.95 10 10 


4.98 10 10 


8.57 10 10 


1.16 10 11 


32 s 


2.2 10 10 


3.69 10 10 


6.25 10 10 


8.4 10 10 


36 Ar 


7.3 10 9 


1.21 10 10 


2.0 10 10 


2.6 10 10 


40 Ca 


1.0 10 10 


1.75 10 10 


2.97 10 10 


3.9 10 10 


44 Ti 


1.3 10 7 


2.02 10 7 


1.0 10 8 


5.2 10 8 


48 Cr 


5.6 10 8 


9.7 10 8 


1.65 10 9 


2.2 10 9 


52 Fe 


2.3 10 9 


5.3 10 9 


1.15 10 10 


1.7 10 10 


56 N j 


1.46 10 10 


4.6 10 10 


1.40 10 11 


2.6 10 11 


60 Zn 


1.0 10 3 


2.7 10 3 


3.27 10 5 


1.5 10 7 
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Table 2: Estimated amounts of IME synthesized during an 
exponential expansion from po = 6.6 x 10 7 g.cirT 3 
down to po = 2 x 10 6 g.cnT 3 with characteristic time t c 







AmiME 




X/M£ 




(s) 


g.cnr 2 


g.cirT 2 


g.crrT 2 






0.05 


9.27 10 10 


6.88 10 10 


1.75 10 10 


0.742 


0.188 


0.1 


1.80 10 11 


1.16 10 11 


5.23 10 10 


0.646 


0.29 


0.2 


3.69 10 11 


1.98 10 11 


1.53 10 11 


0.537 


0.415 


0.3 


5.68 10 11 


2.65 10 11 


2.8 10 11 


0.466 


0.492 
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